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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si)  AND 
METRIC  (SI)  TO  U.  S.  CUSTOMARY  UNITS  OF  MEASUREMENTS 


Units  of  measurement  used  in  this  report  can  be  converted  as 
follows : 


Mult  iply 


JiL 


To  Obtain 


feet 

miles  (U.  S.  statute) 
tons  (2000  lb,  mass) 
miles  (U.  S.  nautical) 


U.  S.  Customary  to  Metric  (SI 

0.30U8 
1.852 
907.I8U7 
1.852 


metres 
kilometres 
kilograms 
kilometres  per  hour 


per  hour 
square  feet  per  second 
feet  per  second  squared 
degrees  (angular) 


0.0929030U     square  metres  per  second 
0.30H8         metres  per  second  squared 
0. 017^5329     radians 


Metres  (Si)  to  U.  S.  Customary 


metres 
kilometres 
square  kilometres 

radians  per  second 

radians 


3.280839       feet 

0.6213711      miles  (U.  S.  statute) 

0.3861021      square  miles  (U.  S. 
statute) 

57.29578        degrees  (angular)  per 
second 

57-29578        degrees  (angular) 


TSUNAMI-WAVE  ELEVATION  FREQUENCY  OF  OCCURRENCE 
IN  THE  HAWAIIAN  ISLANDS 


PART  I:   INTRODUCTION 
Background 

1.  Of  all  water  waves  that  occur  in  nature,  one  of  the  most 
destructive  is  the  tsunami.   The  term  "tsunami,"  originating  from  the 
Japanese  words  "tsu"  (harbor)  and  "nami"  (wave),  is  used  to  describe 
sea  waves  of  seismic  origin.   Tectonic  earthquakes,  i.e.  earthquakes 
that  cause  a  deformation  of  the  sea  bed,  appear  to  be  the  principal 
seismic  mechanism  responsible  for  the  generation  of  tsunamis.   Coastal 
and  submarine  landslides  and  volcanic  eruptions  also  have  triggered 
tsunamis. 

2.  Tsunamis  are  principally  generated  by  undersea  earthquakes  of 
magnitudes  greater  than  6.5  on  the  Richter  scale.   While  the  typical 
height  of  a  tsunami  in  the  deep  ocean  may  be  only  a  foot  or  less,  the 
waves  have  a  tremendous  amount  of  energy  which  is  indicated  by  the  long 
period  between  crests  (5  min  to  several  hours).   Tsunamis  travel  at  the 
shallow-water  wave  celerity  equal  to  the  square  root  of  acceleration  due 
to  gravity  times  water  depth  even  in  the  deepest  oceans  because  of  their 
very  long  wavelengths.   This  speed  of  propagation  can  be  in  excess  of 
500  mph*  in  the  deep  ocean. 

3.  When  tsunami  waves  approach  a  coastal  region  where  the  water 
depth  decreases  rapidly,  wave  refraction,  shoaling,  and  bay  or  harbor 
resonance  may  result  in  significantly  increased  wave  heights.   The  great 
period  and  wavelength  of  tsunami  waves  preclude  their  dissipating  energy 
as  a  breaking  surf;  instead,  they  are  apt  to  appear  as  bores  or  just 
rapidly  rising  water  levels. 

h.      Localized  bathymetric  and  topographic  features  may  cause  the 

*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (Si)  units  and  metric  (Si)  units  to  U.  S.  customary 
units  is  presented  on  page  5- 


effect  of  the  same  tsunami  wave  to  be  drastically  different  from  one 
point  to  another.   For  example,  a  bay  such  as  Hilo  Bay  of  the  island 
of  Hawaii  (Figure  l),  has  a  funneling  effect  that  tends  to  increase  the 
height  of  the  wave.   However,  an  offshore  sandbar  or  reef  such  as  that 
of  Kaneohe  Bay  on  the  island  of  Oahu  (Figure  l)  allows  only  a  diminished 
wave  height  to  reach  the  coast.   It  is  not  impossible  for  the  same  wave 
to  be  50  ft  high  at  one  location  and  only  5  ft  high  a  few  miles  distant 
as  a  result  of  local  factors. 

5.  The  loss  of  life  and  destruction  of  property  due  to  tsunamis 
have  been  immense.   The  Great  Hoei  Tokaido-Nanhaido  tsunami  of  Japan 
killed  30,000  people  in  1707-   In  1868,  the  Great  Peru  tsunami  caused 
25,000  deaths  and  carried  the  frigate  U.S.S.    Waterlee   1,300  ft  inland. 
The  Great  Meiji  Sanriku  tsunami  of  1896  killed  27,122  persons  in  Japan 
and  washed  away  over  10,000  houses.   The  most  recent  major  tsunami  to 
affect  the  United  States,  the  I96U  Alaskan  tsunami,  killed  107  people 
in  Alaska,  k   in  Oregon,  and  11  in  Crescent  City,  California,  and  caused 
over  $100  million  in  damage  on  the  west  coast  of  North  America. 

Tsunamis  in  the  Hawaiian  Islands 

6.  The  Hawaiian  Islands,  a  chain  of  eight  islands  as  shown  in  Fig- 
ure 1,  have  a  history  of  destructive  tsunamis  generated  both  in  distant 
areas  and  locally.   The  earliest  recording  of  a  severe  tsunami  in  the 
Hawaiian  Islands  was  in  1837  when  a  tsunami  from  Chile  reached  an  ele- 
vation of  20  ft  at  Hilo  and  killed  k6   people  in  the  Kau  District  of  the 
island  of  Hawaii  (Figure  l).   Prior  to  l837»  a  number  of  severe  tsunamis 
undoubtedly  reached  the  islands ,  but  unfortunately  no  detailed  records 

were  kept.   Since  l8375  there  have  been  l6  tsunamis  that  have  caused 

2 
significant  damage. 

7.  Most  of  the  destructive  tsunamis  in  the  Hawaiian  Islands  have 
been  generated  along  the  coast  of  South  America,  the  Aleutian  Islands, 
and  the  Kamchatkan  Peninsula.   Approximately  one  fourth  of  all  the 
tsunamis  recorded  in  the  Hawaiian  Islands  have  originated  along  the 
coast  of  South  America,  while  more  than  one  half  have  originated  in  the 
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Kuril-Kamchatka-Aleutian  region  of  the  north  and  northwestern  Pacific. 

8.  Tsunamis  generated  by  local  seismic  events  have  caused  large 
runup  in  the  islands,  especially  on  the  southeast  coast  of  the  big 
island  of  Hawaii.   The  1868  tsunami  produced  the  largest  waves  of  record 
in  the  Hawaiian  Islands  with  60-ft  waves  reported  on  the  South  Puna 
coast  of  the  island  of  Hawaii.   The  most  recent  tsunami  in  Hawaiian 
history  occurred  on  29  November  1975,  when  waves  generated  by  an  earth- 
quake with  an  epicenter  on  the  South  Puna  Coast  reached  elevations  as 
great  as  perhaps  26  ft  along  the  southeast  coast. 

9.  The  most  destructive  tsunami  to  ever  hit  the  islands  in 
terms  of  loss  of  life  and  destruction  of  property  was  the  Great  Aleutian 
tsunami  of  19^-6,  which  killed  173  people  and  produced  waves  over  55  ft 
in  elevation.   Hilo  incurred  $26  million  in  property  damage  attributable 
to  this  tsunami. 

10.  The  i960  Chilean  tsunami  is  the  most  recent  distantly  gen- 
erated tsunami  that  produced  major  effects  in  the  Hawaiian  Islands. 
Sixty-one  lives,  all  at  Hilo,  were  claimed  by  the  tsunami.   Damage 
throughout  the  State  was  estimated  to  be  $23-5  million  of  which  93  per- 
cent occurred  at  Hilo.   Other  major  damage  was  restricted  to  the  Kahului 
area  of  the  island  of  Maui.   Inspection  of  the  damage  at  Hilo  revealed 
much  evidence  of  the  tremendous  forces  developed  by  the  waves.   Twenty- 
ton  boulders  had  been  moved  hundreds  of  feet,  asphaltic  concrete  pave- 
ments were  peeled  from  their  subbase,  and  hundreds  of  automobiles  were 
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tossed  around  and  crushed. 

Purpose  of  Study 

11.  The  purpose  of  this  study  was  to  establish  tsunami-wave  ele- 
vation near  the  shoreline  versus  frequency-of-occurrence  curves  for  the 
Hawaiian  Islands  and  to  recommend  future  efforts  to  determine  methods  of 
calculating  runup  for  those  areas  where  substantial  flooding  of  low- 
lying  areas  causes  runup  not  to  be  equal  to  wave  elevation  near  the 
shoreline.   This  information  is  required  by  the  U.  S.  Army  Engineer 
Division,  Pacific  Ocean,  for  use  in  tsunami  flood  hazard  evaluations 


for  flood/plain  management  and  flood  insurance  rate  calculations.   The 
results  of  this  study  should  not  be  used  to  determine  evacuation  zones 
based  upon   l-in-100-year  inundation   levels  without  adding  a  safety 
factor  through  risk  oalcualtions    (Part  \V) .    The  odds  that  an  elevation 
greater  than  the  general  l-in-100-year  elevation  will  occur  within  a 
short  time  span  are  not  negligible.   Any  land  development  that  would 
expose  human  life  to  possible  danger  should  be  initiated  only  after  an 
evaluation  of  the  possible  risk. 
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PART  II:   APPROACH 
Time  Period  Analyzed 

12.  It  is  necessary  to  use  historical  data  of  tsunami  occurrence 
covering  the  greatest  possible  time  span  to  properly  determine  frequency 
of  occurrence  of  tsunami  elevations  for  the  Hawaiian  Islands.   Tsunami 
activity  has  not  been  uniform  in  these  islands  during  recorded  history. 
For  example,  the  two  largest  and  four  of  the  ten  largest  tsunamis  strik- 
ing Hilo  from  1837  through  1976  occurred  during  the  15-year  period  from 
19^6  through  i960.   Two  of  the  tsunamis  from  19^6  through  i960  originated 
in  the  Aleutian  Islands,  one  in  Kamchatka,  and  one  in  Chile.   However, 
six  of  the  ten  largest  tsunamis  occurred  during  the  109-year  period 

from  1837  through  19^5  with  three  originating  in  Chile,  two  in  Kanchatka, 
and  one  in  Hawaii.   Therefore,  both  the  frequency  of  occurrence  and 
place  of  origin  of  tsunamis  have  been  remarkably  variable.   Any  study, 
such  as  Reference  h,   using  a  short  time  span  that  includes  the  period 
from  19^+6  through  i960  will  predict  a  significantly  more  frequent  oc- 
currence of  large  tsunamis  than  is  warranted  by  historical  data  from 
1837  through  1976.   Although  the  quantitative  accuracy  of  the  data  for 
tsunamis  from  1837  through  19^5  may  be  questionable,  there  is  little 
doubt  that  the  recorded  occurrence  of  large  tsunamis  is  accurate  (i.e., 
tsunamis  noted  as  being  significant  were  indeed  so,  and  major  tsunamis 
did  not  occur  and  go  unrecorded). 

13.  From  an  analysis  of  tsunami  data  for  Hilo  the  errors  intro- 
duced in  frequency-of-occurrence  calculations  by  consideration  only 

of  a  short  period  that  includes  the  unrepresentative  years  from  19^6 
through  i960  will  be  greater  than  the  errors  resulting  from  possible 
observational  inaccuracies  of  the  nineteenth  century.   A  calculated 
l-in-100-year  tsunami  elevation  for  Hilo  is  27.3  ft.   This  was  deter- 
mined by  a  least-squares  fit  of  the  data  using  the  logarithmic  fre- 
quency distribution  discussed  later  and  based  upon  data  compiled  by 
Cox  for  the  ten  largest  tsunamis  in  Hilo  from  1837  through  1976.   (His 
data  were  actually  through  196U,  but  no  large  events  have  occurred  since 
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I96U) .   Increasing  by  50  percent  the  elevations  for  five  of  the  ten  large 
tsunamis  occurring  during  the  nineteenth  century  yields  a  l-in-100-year 
elevation  of  30. h   ft.   The  l-in-100-year  elevation  that  is  based  just 
upon  the  large  tsunamis  during  the  period  of  accurate  survey  measurements 
in  Hilo  from  19^6  through  1976  is  UU.2  ft.   Since  the  largest  elevation 
in  Cox's  data  for  the  1^0-year  period  from  183T  through  1976  was  28  ft 
for  the  i960  Chilean  tsunami,  the  HU.2-ft  elevation  is  obviously  much 
too  large. 

ik.      Since  an  overestimate  of  the  frequency  of  major  tsunamis  at 
Hilo,  which  has  the  most  complete  data  for  tsunamis  in  the  Hawaiian 
Islands,  did  result  from  using  an  analysis  based  upon  the  short  period 
from  19^-6  through  I960,  it  was  concluded  that  overestimates  will  occur 
at  all  other  locations  throughout  the  islands  from  analyses  based  upon 
similar  short  periods.   The  exceptionally  frequent  occurrence  of  major 
tsunamis  in  Hilo  from  19^6  through  i960  is  a  property  of  the  unusual 
activity  of  tsunami  generation  areas  and  not  of  special  properties 
of  Hilo. 

15.  The  unusual  tsunamigenic  activity  in  generation  regions  is 
reflected  in  the  most  up-to-date  catalog  of  tsunami  occurrence  in  the 
circumpacific  area  compiled  by  Soloviev.    This  catalog  shows,  for 
example,  that  since  1788  the  three  tsunamis  generated  with  the  greatest 
intensity  in  the  Aleutian- Alaskan  region  occurred  in  19^-6,  1957,  and 
196U.   No  significant  tsunamis  were  generated  during  the  period  from 
1837  through  19^-5  in  this  region  and  there  also  were  no  reports  of 

significant  tsunamis  in  the  Hawaiian  Islands  which  might  have  originated 

2 
in  this  region.   However,  Pararas-Carayannis  reports  a  minor  tsunami 

arriving  in  the  Hawaiian  Islands  from  the  Aleutian-Alaskan  region  on 

10  November  1938  that  produced  an  elevation  in  Hilo  of  0.3  m  and  in 

Honolulu  of  0.1  m.   This  tsunami  also  is  reported  as  a  minor  one  in 

Soloviev' s  catalog.   Clearly,  tsunami  generation  was  much  more  active 

in  the  Aleutian-Alaskan  region  from  19^+6  through  196U  than  from  1837 

through  19^5. 

16.  The  19^6  and  1957  tsunamis  are  the  dominant  tsunamis  for 
northern  coasts  of  many  of  the  Hawaiian  Islands  (the  19^+6  tsunami  is 
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one  of  the  largest  tsunamis  for  most  locations  even  on  southern  coasts). 
For  example,  at  Haena,  Kauai,  the  1957  tsunami  elevation  was  52.5  ft; 
the  19^+6  tsunami  elevation,  45.0  ft;  and  the  i960  Chilean  tsunami  eleva- 
tion, 9-2  ft.   A  least-squares  fit  of  these  data  using  the  logarithmic 
distribution  discussed  later  for  the  period  19^6  through  1976  yields  a 
l-in-100-year  elevation  of  97.6  ft.   However,  since  Aleutian- Alaskan 
tsunamis  clearly  produce  the  largest  elevations  at  Haena  (and  on  all  of 
the  northern  Kauai  coast),  the  52.5-ft  elevation  must  be  the  largest 
tsunami  elevation  at  Haena  since  at  least  1788.   As  reported  by  Solo- 
viev,   the  1788  tsunami  generated  in  the  Aleutian- Alaskan  region  was 
quite  large  but  produced  unknown  consequences  to  the  Hawaiian  Islands. 
A  logarithmic  distribution  fit  of  the  same  data  for  the  years  1837 
through  1976  yields  a  l-in-100- year  elevation  of  hk.5   ft.   Therefore, 
the  period  from  19^6  through  1976  is  clearly  unrepresentative  of  long- 
period  tsunami  activity  in  the  Hawaiian  Islands  due  to  Aleutian-Alaskan 
generated  tsunamis. 

17-   Soloviev's  catalog  also  lists  the  i960  Chilean  tsunami  as 
having  the  largest  intensity  of  any  tsunami  generated  in  South  America 
in  recorded  history  (4l5-year  period  from  1562  through  1976).   Reporting 
of  tsunamis  on  the  west  coast  of  South  America  is  quite  good  with  105 
tsunamis  reported  during  the  415-year  period.   A  recent  revision  of  the 
Richter  scale  for  measuring  earthquake  energy  by  H.  Kanamori  of  the 
California  Institute  of  Technology*  and  a  consequent  revision  of  Richter 
magnitude  estimates  for  several  important  large  earthquakes  listed  a  new 

magnitude  of  9.5  for  the  i960  Chilean  earthquake,  thus  making  this 

7 

earthquake  the  largest  in  recorded  history.   McGarr  recently  calculated 

seismic  moments  of  earthquakes  and  found  that  earthquakes  much  larger 
than  the  Chilean  earthquake  of  i960  do  not  seem  possible  in  view  of  the 
bounds  associated  with  the  relative  motion  of  tectonic  plates.   Earth- 
quakes as  large  as  the  i960  event  are  probably  quite  rare.   Plafker  and 

o 

Rubin,   for  example,  claim  that  relative  submergence  occurred  along  the 
Alaskan  coast  for  at  least  900  years  before  the  1964  Alaskan  earthquake. 

*  Associated  Press  news  article,  March  1977. 
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The  period  from  19^6  through  i960  certainly  cannot  be  considered  as  being 
representative  of  long-period  tsunami  activity  in  South  America  when  the 
greatest  tsunami  (and  earthquake)  in  recorded  history  for  this  region 
occurred  in  i960. 

18.  Tsunamis  arriving  in  the  Hawaiian  Islands  from  Kamchatka  may 
not  be  too  abnormally  concentrated  within  the  period  from  19^+6  through 
i960.   However,  tsunamis  from  Kamchatka  rarely  produce  the  dominant 
flood  elevations  for  locations  in  the  Hawaiian  Islands.   Two  large 
tsunamis  from  Kamchatka  arrived  in  the  Hawaiian  Islands  during  the  years 
1837  through  19^5  (return  period  approximately  55  years)  and  one  during 
the  years  19^6  through  1976  (return  period  approximately  31  years). 
This  difference  in  return  periods  may  not  be  significant,  but  again  it 
points  to  the  19^-6  through  i960  period  as  not  being  representative  of 
long-term  activity.   It  is  interesting  to  note  that  Soloviev  reports 
two  large  tsunamis  generated  in  Kamchatka  from  1737  through  1836  (return 
period  approximately  51  years).   Thus,  the  return  period  for  large 
tsunamis  was  approximately  the  same  for  the  periods  1737  through  1836 
and  1837  through  19^5. 

19.  Evidence  that  tsunamis  in  the  Hawaiian  Islands  were  unusually 
active  from  19^-6  through  i960  also  exists  at  Kahului  and  Lahaina,  Maui. 
Kahului,  Lahaina,  Hilo,  and  Honolulu  are  the  only  locations  in  the 
islands  with  historical  data  as  far  back  as  1837-   Four  of  the  five 
largest  tsunami  elevations  recorded  at  Kahului  occurred  during  the 
years  19^6  through  i960.   Three  of  the  four  largest  tsunami  elevations 
recorded  at  Lahaina  occurred  during  the  same  period.   In  addition,  the 
19^-6,  1957,  and  i960  tsunamis  produced  the  three  largest  elevations 
recorded  anywhere  in  the  Hawaiian  Islands  since  1837. 

20.  The  main  evidence  that  the  period  19^6  through  1976  might 

not  be  unusual  for  locations  other  than  Hilo  is  the  tsunami  data  at 

Honolulu.   The  greatest  elevation  in  Honolulu  used  in  Reference  k   is  a 

2.U-m  elevation  for  the  1837  tsunami  from  Chile.   However,  this  2.k   m  is 

not  an  elevation  but  a  drop  in  water  level.   According  to  Pararas- 

2 
Carayannis,   "at  5:00  the  water  began  receding  to  2.4  m  leaving  Honolulu 

Harbor  partly  dry,  then  slowly  returned.   Wave  action  lasted  until  the 
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following  day."   There  is  no  indication  of  the  maximum  crest  elevation 
achieved  other  than  the  fact  that  it  was  less  than  2.U  m.   The  amplitude 
of  troughs  is  of  no  importance  to  study  of  elevation  frequencies  and  is 
not  comparable  to  data  just  composed  of  crest  amplitudes.   Indeed,  a 
trough  of  the  i960  tsunami  from  Chile  measured  in  Honolulu  may  have  had 
an  amplitude  as  great  as  twice  the  maximum  crest  amplitude  (Figure  9, 
page  UO),  yet  Reference  k   just  used  the  reported  crest  amplitude  for 
this  tsunami.   If  the  1837  tsunami  is  neglected,  the  largest  elevation 
produced  in  Honolulu  by  distantly  generated  tsunamis  occurred  during 
the  19^6  Aleutian  and  the  i960  Chilean  tsunami.   Again,  events  during 
the  years  19^+6  through  i960  dominate  the  maximum  flood  elevations. 

21.  Tsunami  data  for  Honolulu  must  be  carefully  considered  before 
wide-ranging  conclusions  based  upon  the  data  are  made.   Since  tsunami 
elevations  are  small  at  Honolulu,  the  tide  level  may  have  a  significant 
influence  on  maximum  elevations.   Thus,  it  is  difficult  to  compare 
probable  tsunami  crest  amplitudes  for  historical  tsunamis  because  the 
reported  elevations  may  have  significant  tidal  components.   Also,  the 
elevation  used  in  Reference  k   for  the  l8Ul  tsunami  is  a  trough  amplitude, 
and  elevations  for  the  1868  and  1877  tsunamis  from  Chile  are  not  known 
for  Honolulu. 

22.  Recently,  D.  C.  Cox*  has  suggested  the  possibility  of  tsunamis 
generated  in  the  central  or  western  Aleutians  having  occurred  during  the 
period  from  l83T  through  19^5  and  produced  significant  elevations  in  the 
Hawaiian  Islands  at  locations  other  than  the  populated  communities  of 
Hilo,  Honolulu,  Kahului,  and  Lahaina.   As  evidence,  Cox  cites  the  wave 
damage  at  Haleiwa,  Oahu,  in  1872  probably  resulting  from  a  10- ft  eleva- 
tion of  water.   However,  he  further  states  that  the  damage  might  have 
been  caused  by  a  tsunami  from  the  north,  but  it  was  more  likely  a  re- 
sult of  storm  waves.   A  tsunami  later  in  1872,  possibly  from  the  north, 
resulted  in  rises  of  k   ft  at  Hilo,  1  ft  at  Honolulu,  1.5  ft  at  Nawiliwili, 
Kauai,  and  probably  about  10  ft  at  Hanalei,  Kauai.   An  1878  tsunami, 

also  possibly  from  the  north,  resulted  in  a  reported  runup  of  10  ft  at 

*  Personal  letter,  7  March  1977. 
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Waialua,  Oahu,  and  probably  about  10  or  15  ft  at  Maliko,  Honomanu,  and 
Halehaku  on  the  island  of  Maui.   Cox  believes  that  these  events  suggest 
that  Aleutian-Alaskan  tsunamis  might  have  been  more  active  at  some  loca- 
tions in  the  Hawaiian  Islands  from  1837  through  19^+5  than  indicated  in 
this  study. 

23.   The  evidence  that  no  major  tsunamis  originating  in  the 
Aleutian -Alaskan  area  from  1837  through  19^5  occurred  and  went  unreported 
is  persuasive.   It  is  true  that  Hilo  is  not  as  strongly  responsive  to 
tsunamis  generated  in  the  western  Aleutians  as  it  is  to  tsunamis  gen- 
erated in  the  central  Aleutians  or  the  eastern  Alaskan  area.   However, 
the  1957  tsunami  from  the  western  Aleutians,  although  not  large  compared 
with  many  of  the  historical  tsunamis  in  Hilo,  did  produce  an  elevation 

only  1  ft  less  than  that  of  the  1868  Chilean  tsunami,  which  was  char- 

2 
acterized  as  having  caused  "severe  damage"  in  Hilo.    During  the  1957 

2 
tsunami,  "buildings  along  the  waterfront  were  badly  damaged;"   Coconut 

Island  in  Hilo  Harbor  was  covered  by  3  ft  of  water  and  the  bridge  to  the 
island  was  destroyed.   Kahului ,  Maui,  also  is  sensitive  to  Aleutian- 
Alaskan  tsunamis  (since  it  faces  north,  as  does  Hilo)  and  has  historical 
data  back  to  1837-   The  19^6  Aleutian  tsunami  was  very  destructive  in 
Kahului ,  and  even  the  1957  tsunami  produced  an  elevation  the  same  as 
that  of  the  1868  Chilean  tsunami,  which  caused  considerable  damage  in 
Kahului.   Yet,  like  Hilo,  there  is  no  report  of  a  significant  Aleutian- 
Alaskan  tsunami  during  the  period  1837  through  19^5-   Furthermore, 
Honolulu  and  Lahaina  were  sufficiently  active  ports  in  the  nineteenth 
century  to  notice  water  level  changes  of  only  1  m  during  the  18^-1  Kam- 
chatkan  tsunami  (one  of  the  smallest  major  tsunamis  in  the  Hawaiian  Is- 
lands during  the  nineteenth  century).   Again,  no  significant  Aleutian- 
Alaskan  tsunamis  during  the  years  1837  through  19^+5  were  reported  at 
these  locations. 

2h.      There  also  is  evidence  from  the  Aleutian-Alaskan  area  that 
the  possible  events  of  1872  and  1878  were  not  tsunamis  generated  in  this 
area.   George  Davidson,  U.  S.  Coast  and  Geodetic  Survey,  was  sent  to 

Q 

Alaska  in  the  summer  of  1867  to  lead  a  survey  team.    The  party  estab- 
lished a  tide  gage  at  Kodiak,  which  was  sensitive  enough  to  record  the 
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1868  tsunami  from  Chile  and  the  1883  blast  waves  from  Krakatoa,  Indo- 
nesia.  (The  1883  event  was  recorded  by  a  tide  gage  in  Honolulu  as 
having  a  height  of  0.2  m. )  Apparently,  Davidson  made  no  report  of 
tsunamis  generated  in  1872  or  1878  in  the  Alaskan  area. 

25.  The  elevations  reported  by  Cox  for  the  possible  tsunamis  of 
1872  and  1878  are  not  particularly  large  compared  with  other  historical 
events  at  the  indicated  locations.   For  example,  at  least  four  tsunamis 
at  Haleiwa,  Oahu,  just  during  the  19-year  period  from  19^+6  through  196^ 
had  elevations  greater  than  the  possible  tsunami  of  1872.   Actually,  if 
the  possible  tsunamis  of  1872  and  1878  are  added  to  the  analysis  de- 
scribed later  in  this  report  at  the  indicated  locations,  the  l-in-100- 
year  elevations  are  only  changed  by  fractions  of  a  foot.   Since  these 
possible  events  do  not  significantly  influence  the  l-in-100-year  eleva- 
tions for  locations  on  Kauai,  Oahu,  Maui,  or  for  Hilo,  Hawaii,  they  are 
probably  not  significant  for  elevation  calculations  elsewhere  in  the 
Hawaiian  Islands.   Furthermore,  the  possible  occurrence  of  these  events 
is  not  an  indication  that  the  Aleutian-Alaskan  tsumanis  were  active 
enough  at  locations  in  the  Hawaiian  Islands  during  the  years  1837  through 
19^+5  for  the  period  from  19^+6  through  i960  to  be  representative  of  long- 
period  Aleutian-Alaskan  tsunami  activity. 

Interpolation  of  Historical  Data 

26.  Historical  data  of  tsunami  activity  in  the  Hawaiian  Islands 
are,  of  course,  often  limited  to  certain  locations.   Information  on 
tsunami  activity  in  the  islands  prior  to  the  19^-6  tsunami  is  concen- 
trated in  Hilo,  Hawaii,  and  Honolulu,  Oahu,  and  to  a  lesser  extent  in 
cities  such  as  Kahului  and  Lahaina  on  the  island  of  Maui.   Even  data  for 
tsunamis  from  19^+6  through  1976  are  absent  or  fragmentary  (i.e.,  data 
exist  only  for  certain  of  the  events)  for  much  of  the  coastline  of  the 
islands.   Therefore,  it  is  necessary  to  rely  on  more  than  just  available 
historical  data  to  determine  tsunami  occurrence  frequencies.   In  this 
study,  a  hybrid  finite  element  numerical  model  was  used  to  supplement 
historical  information  by  allowing  for  interpolation  between  historical 
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data.   The  model  is  described  in  Part  III  and  verified  by  a  comparison 
of  the  numerical  model  calculations  and  historical  data  recorded  during 
the  i960  and  196^  tsunamis. 

27.  The  numerical  model  is  used  to  fill  in  historical  data  gaps 
for  tsunamis  from  19^-6  through  196U  by  providing  relative  responses  of 
the  Hawaiian  Islands  to  tsunamis.   Although  the  deepwater  wave  form  of 

a  tsunami  such  as  that  of  19^6,  for  example,  is  not  known,  the  direction 
of  approach  of  this  tsunami  and  its  range  of  wave  periods  are  known. 
By  inputting  sinusoidal  waves  of  unit  amplitude  from  the  same  direction 
as  the  19^6  tsunami  into  the  numerical  model  over  a  band  of  wave  periods, 
the  average  relative  response  of  the  islands  to  a  tsunami  similar  to  the 
19^6  tsunami  is  determined.   If  historical  data  exist  at  one  location, 
wave  elevations  at  a  nearby  location  for  which  historical  data  do  not 
exist  are  determined  by  multiplying  the  historical  data  at  the  first 
location  by  the  ratio  of  the  response  calculated  by  the  numerical  model 
at  the  second  location  and  the  response  calculated  at  the  first.   The 
numerical  model  takes  into  account  the  major  processes  that  would  cause 
different  wave  elevations  at  the  two  locations.   That  is,  the  model 
calculates  shoaling,  refraction,  diffraction,  reflection,  resonance, 
shielding  of  the  back  side  of  an  island  by  the  front  side,  and  inter- 
actions between  islands.   Historical  data  from  19^-6  through  196^  were 
taken  from  the  most  recent  compilation  of  these  data.    Part  III  de- 
scribes the  use  of  the  numerical  model  for  interpolation  in  detail. 

28.  As  an  example  of  the  use  of  the  numerical  model  to  fill 
historical  data  gaps  for  tsunamis  from  19^-6  through  196U,  consider  the 
19^6  tsunami  that  was  recorded  on  the  island  of  Lanai  only  at  Kaumalapua 
Harbor  and  Manele  Bay.   Wave  elevations  can  be  calculated  for  the  19^-6 
tsunami  at  other  locations  on  the  coastline  of  Lanai  by  multiplying 
these  historical  values  by  the  ratio  of  the  response  calculated  at  lo- 
cations on -the  coastline  by  the  numerical  model  and  the  response  calcu- 
lated at  Kaumalapua  Harbor  and  Manele  Bay  by  the  numerical  model.   Thus, 
wave  elevations  for  the  19^6  tsunami  can  be  calculated  for  all  of  Lanai. 

29.  Some  islands  lack  any  historical  data  on  a  particular  tsu- 
nami.  For  example,  there  are  no  data  for  the  i960  tsunami  on  Lanai. 
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Historical  data  for  such  cases  can  be  reconstructed  by  using  historical 
data  from  areas  of  nearby  islands  in  conjunction  "with  the  numerical  model. 
For  example,  the  west  coast  of  Maui  (near  Lahaina)  is  less  than  10  miles 
from  the  east  coast  of  Lanai  and  has  historical  data  for  the  i960 
tsunami.   Wave  elevations  can  be  calculated  on  Lanai  for  the  i960  tsunami 
by  multiplying  the  historical  wave  elevations  on  the  west  coast  of  Maui 
by  the  ratio  of  the  response  calculated  on  Lanai  by  the  numerical  model 
and  the  response  calculated  on  the  west  coast  of  Maui  by  the  numerical 
model.   The  wave  elevations  calculated  for  the  i960  tsunami  on  Lanai 
using  this  approach  are  consistent  with  the  qualitative  reports  on 
tsunami  activity  during  the  i960  tsunami.   Damages  on  Lanai  during  the 
i960  tsunami  were  slight  ($2000)  and  confined  to  Kaumalapau  Harbor  where 
three  small  boats  were  sunk.   One  beach  house  was  damaged  on  the  island. 
The  reconstructed  wave  elevations  for  the  i960  tsunami  on  Lanai  were 
fairly  small.   They  were  smaller  than  the  elevations  recorded  on  the 
west  coast  of  Maui  where  the  i960  tsunami  caused  damage  at  Lahaina 
amounting  to  $17,000. 

Reconstruction  of  Historical  Data 

30.   Tsunamis  originating  near  the  Aleutian  Islands,  Kamchatka, 
and  Chile  occurred  from  19^+6  to  196U.   Therefore,  the  response  is  known 
of  many  areas  in  the  Hawaiian  Islands  to  tsunamis  originating  in  the 
three  main  locations  where  tsunamis  of  destructive  power  in  these 
islands  have  historically  been  generated.   It  is  assumed  in  this  study 
that  tsunamis  generated  in  a  single  source  region  (Kamchatka  or  Chile 
but  not  the  Aleutians)  approach  the  islands  from  approximately  the  same 
direction  and  have  energy  lying  in  the  same  band  of  wave  periods.   The 
difference  in  wave  elevations  at  the  shoreline  in  the  Hawaiian  Islands 
produced  by  tsunamis  generated  at  different  times  in  the  same  region  is 
attributed  mainly  to  differences  in  deepwater  wave  amplitudes.   For  ex- 
ample, the  18^1  tsunami  from  Kamchatka  produced  a  wave  elevation  in  Hilo 

that  was  approximately  25  percent  greater  than  that  of  the  1952  tsunami 

2 
from  Kamchatka.    The  same  relative  magnitude  of  the  two  tsunamis  is 

used  for  all  of  the  islands  to  determine  the  elevation  that  must  have 
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occurred  in  l8Ul  at  some  location,  knowing  the  elevation  which  did  occur 
in  1952.   Therefore,  knowing  the  elevations  of  tsunamis  at  a  location 
from  19U6  to  i960  and  the  response  of  Hilo  to  tsunamis  from  1837  to  i960 
allows  a  reconstruction  of  the  elevations  of  tsunamis  that  occurred 
prior  to  19^6  but  were  not  recorded. 

31.  The  assumption  that  tsunamis  generated  in  Kamchatka  and  Chile 
approach  the  Hawaiian  Islands  from  nearly  the  same  direction  is  justi- 
fied by  the  small  spatial  extent  of  the  known  generation  areas  in  Kam- 
chatka and  results  of  a  study   that  considered  the  propagation  of 
tsunamis  from  Chile  using  a  finite  difference  numerical  model.   This 
study  indicated  that  the  directional  effects  for  tsunamis  originating 
along  the  Chilean  coast  are  small  in  the  Hawaiian  Islands  (probably 
because  the  generation  areas  in  Chile  subtend  a  relatively  small  angle 
with  respect  to  the  islands).   The  position  of  the  Aleutian-Alaskan 
Trench  relative  to  the  Hawaiian  Islands  does  introduce  important  direc- 
tional efforts  for  tsunamis  generated  in  the  Aleutian-Alaskan  area. 
However,  these  effects  are  known  from  historical  observations  for  tsu- 
namis generated  in  the  western  Aleutians  (1957  tsunami),  central 
Aleutians  (19^6),  and  eastern  Alaskan  area  (196U). 

32.  Many  observers  have  noted  that  wave  elevations  produced  at  a 
location  by  "tsunamis  of  diverse  geographic  origin  are  strikingly  dif- 
ferent, whereas  those  from  nearly  the  same  origin  are  remarkably 
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similar."    A  recent  study   shows  that  the  normalized  frequency  dis- 
tribution for  tsunamis  of  the  same  origin  have  the  same  distribution. 
Therefore,  the  assumption  is  quite  reasonable  that  heights  of  past 
tsunamis  for  which  no  data  exist  can  be  estimated  by  considering  tsu- 
namis of  the  same  origin  for  which  data  do  exist. 

33.  In  his  report,  Cox  noted  the  problem  of  selecting  elevations 
measured  during  surveys  of  runup  and  inundation  for  tsunamis  from  19^6 
through  I96U  in  Hilo  Bay  that  are  comparable  with  the  earlier  visual 
elevations  recorded  there  from  1837  to  19^6.   He  claimed  that  it  could 
be  assumed  that  the  reporters  of  early  tsunamis  did  not  record  the  mini- 
mum level  reached,  nor  even  an  average.   Furthermore,  he  believed  that 
it  could  not  be  assumed  that  they  identified  the  absolute  maximum  in  the 
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absence  of  surveys  comparable  with  the  ones  from  1946  through  1964. 
Therefore,  in  order  to  make  the  recorded  values  from  1946  through  1964 
comparable  with  earlier  historical  tsunamis,  he  chose  as  the  "normal 
maximum"  elevation  the  upper  decile  value  for  each  of  the  tsunamis  from 
1946  through  1964  from  inundation  heights  measured  along  the  shores  of 
Hilo  Bay  between  the  west  bank  of  the  Wailuku  River  and  the  root  of  the 
breakwater.   His  compilation  of  historical  data  and  comparable  adjusted 
values  were  used  to  determine  the  ratios  needed  to  reconstruct  the  ele- 
vation of  tsunamis  that  occurred  prior  to  1946  at  locations  lacking  such 
historical  data. 

34.   Tsunami  elevations  along  the  coasts  of  the  Hawaiian  Islands 
can  be  calculated  for  the  ten  largest  tsunamis  in  Hilo  from  1837  through 
1976  by  using  the  finite  element  numerical  model  to  fill  in  data  gaps 
for  those  tsunamis  of  the  ten  that  occurred  from  1946  through  1964  and 
then  taking  these  results  and  the  ratios  from  Cox's  study  to  reconstruct 
the  elevations  of  those  tsunamis  that  occurred  prior  to  1946.   For  much 
of  the  coast  of  the  Hawaiian  Islands,  the  ten  largest  tsunamis  since 
1837  would  be  the  same  as  the  ten  largest  tsunamis  in  Hilo.   That  is, 
the  ten  largest  would  be  the  i960,  1946,  1923,  1837,  1877,  l84l,  1957, 
1952,  and  the  two  in  1868.   Of  course,  the  order  of  the  ten  largest 
would  vary  from  location  to  location.   For  example,  the  largest  tsunami 
on  the  South  Puna  coast  of  the  island  of  Hawaii  was  the  locally  gen- 
erated tsunami  of  1868.   There  are  some  locations  where  one  or  more  of 
the  ten  largest  tsunamis  to  have  occurred  are  not  among  the  ten  largest 
in  Hilo  from  1837  through  1976.   For  example,  the  1896  tsunami  generated 
near  Japan  is  the  largest  to  hit  Keauhau  on  the  Kona  coast  of  the  big 
island  of  Hawaii.   The  ten  largest  tsunamis  at  each  location  are  de- 
termined by  adding  such  historical  tsunamis  and  eliminating  a  corres- 
ponding number  of  the  smallest  tsunamis  from  the  initial  compilation 
at  each  location  of  the  ten  that  were  the  largest  at  Hilo.   The  follow- 
ing tabulation  lists  the  sixteen  tsunamis  considered  in  determining 
the  ten  largest  tsunamis  for  locations  along  the  Hawaiian  coastline, 
the  source  location  of  the  tsunami,  and  the  islands  experiencing 
significant  wave  elevations: 
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Island(s) 

Data 

Source 

Affected 

7  November  1837 

Chile 

All 

IT  May  l8Ul 

Kamchatka 

All 

2  April  1868 

Hawaii 

Hawaii 

13  August  1868 

Chile 

All 

10  May  l8TT 

Chile 

All 

15  June  I896 

Japan 

Hawai  i 

16  August  1906 

Chile 

Maui 

30  April  1919 

Tonga- Island 

Hawai  i 

3  February  1923 

Kamchatka 

All 

2  March  1933 

Japan 

Hawaii 

1  April  19^6 

Aleutian  Islands 

All 

k   November  1952 

Kamchatka 

All 

9  March  1957 

Aleutian  Islands 

All 

23  May  i960 

Chile 

All 

27  March  196k 

Alaska 

All 

29  November  1975 

Hawaii 

Hawai  i 

35.  Historical  data  for  the  tsunamis  prior  to  19^+6  listed  above 
were  taken  from  Reference  2.   Tsunamis  listed  as  being  questionable  in 
Reference  2  are  not  listed.   Data  for  the  29  November  1975  tsunami  were 
taken  from  Reference  ik.      These  data  were  used  in  conjunction  with  data 
from  Reference  2  on  the  2  April  1868  tsunami  to  fill  in  data  gaps  for 
the  1868  tsunami.   Both  the  1868  and  1975  tsunamis  were  generated  on  the 
southeastern  coast  of  the  island  of  Hawaii.   Wave  elevations  in  the  Kau 
District  at  cities  such  as  Honuapo  and  Punaluu  were  similar  for  the  two 
tsunamis.   Wave  elevations  along  the  South  Puna  coast,  however,  were 
much  larger  for  the  1868  tsunami  than  the  1975  tsunami.   The  1868  wave 
elevations  were  reconstructed  for  the  South  Puna  coast  by  multiplying  the 
1975  elevations  by  ratios  determined  at  locations  such  as  Keauhau  Land- 
ing on  the  South  Puna  coast  where  data  exist  for  both  the  1868  and  1975 

t  sunami  s . 

Frequency  of  Occurrence  Distribution 

36.  Cox  found  that  the  logarithm  of  tsunami  frequency  of  occur- 
rence is  linearly  related  to  tsunami  elevations  for  the  ten  largest 
tsunamis  occurring  from  1837  to  I96U  in  Hilo.    This  result  is  not 
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unexpected  considering  the  fact  that  Soloviev  of  the  Soviet  Union  has 
shown  a  similar  relationship  between  tsunami  frequency  of  generation 
and  the  intensity  of  the  tsunamis  near  the  generation  area  for  moderate 
to  large  tsunamis.    Also,  earthquake  intensity  and  frequency  of  occur- 
rence have  been  similarly  related  by  Gutenberg  and  Richter.    Further- 

17 
more,  Wiegel   found  the  same  relationship  for  historical  tsunamis  at 

Hilo,  Hawaii;  San  Francisco,  California;  and  Crescent  City,  California; 

-i  o 

and  Adams,    for  tsunamis  at  Kahuku  Point,  Oahu.   A  recent  study  by 

19 
Rascon  and  Villarreal   revealed  this  same  relationship  for  historical 

tsunamis  on  the  west  coast  of  Mexico  (data  from  1732)  and  on  the 
Pacific  West  coast  of  America,  excluding  Mexico. 

37.   After  the  ten  largest  tsunamis  from  1837  through  1976  are  de- 
termined at  locations  all  along  the  coasts  of  the  Hawaiian  Islands,  el- 
evation versus  frequency  of  occurrence  is  fit  at  each  location  using  the 
least-squares  techniques  by  curves  that  can  be  represented  by  the 
equations 


h  =  -B  -  A  log1QF  (1) 


where 


h*  =  elevation  of  maximum  tsunami -wave  crest  above  mean  sea  level 
200  ft  shoreward  of  the  coastline 

F  =  frequency  per  year  of  occurrence 

20 
According  to  Cox,   unless  there  is  evidence  to  the  contrary,  recorded 

historical  runup  heights  may  reasonably  be  considered  to  have  been  mea- 
sured 200  ft  inland  of  the  coastline.   Since  the  coefficients  of  Equa- 
tion 1  are  based  ultimately  on  recorded  historical  data,  the  values  of 
h  determined  from  this  equation  are  elevations  200  ft  shoreward  of  the 
coastline.   He  fit  his  Hilo  data  by  a  curve  that  can  be  represented  by 
a  similar  equation.   Plots  of  the  coefficients  A  and  B  versus  loca- 
tion along  the  coasts  of  all  the  islands  are  shown  in  Part  IV.   The 
negative  signs  in  Equation  1  were  used  so  that  these  coefficients 
would  be  positive.   Furthermore,  these  coefficients  can  be  found  for  any 
locality  of  the  Hawaiian  Islands  and  used  in  Equation  1  to  calculate 

*  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  de- 
fined in  the  Notation  (Appendix  A). 
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"h"  for  any  desired  "F"  within  a  certain  range  of  frequencies. 

38.   Cox  found  that  the  linear  relation  between  the  logarithm  of 
tsunami  frequency  of  occurrence  and  tsunami  elevations  at  Hilo  was  valid 
for  frequencies  of  occurrence  as  high  as  0.1  per  year  (l-in-10-year 
tsunami).   He  found  that  the  frequency-elevation  relation  for  higher 
frequencies  followed  approximately  a  power  law.   Equation  1  should, 
therefore,  he  used  only  for  frequencies  of  occurrence  lower  than  0.1  per 
year.   There  are  some  locations,  however,  where  the  linear  relation  is 
valid  only  for  frequencies  of  occurrence  lower  than  0.05  (l-in-20-year 
tsunami)  per  year.   This  problem  is  discussed  in  Part  IV,  and  a  table 
(Table  l)  is  presented  that  gives  0.1-per-year  elevations  for  the  entire 
coastline  of  the  Hawaiian  Islands.   It  is  recommended  that  Table  1  be 
used  for  the  0.1-per-year  elevations  and  Equation  1  for  frequencies 
lower  than  or  equal  to  0.05  per  year. 

39 •  Equation  1  was  used  as  a  frequency  of  occurrence  distribution 
because  it  has  been  found  to  agree  with  historical  data  at  several  loca- 
tions, as  discussed  earlier.  Other  distributions  may  agree  with  histor- 
ical data  equally  well.   For  example,  the  Gumbel  distribution  has  been 

21 
used  in  the  past  to  study  annual  streamflow  extremes.    Borgman  and 

22 
Resio   illustrate  the  use  of  this  distribution  to  determine  frequency 

curves  for  nonannual  events  in  wave  climatology.   To  investigate  the 
sensitivity  of  calculations  of  l-in-100-year  elevations  on  the  assumed 
frequency  distribution,  the  approach  of  Borgman  and  Resio  was  applied  to 
the  Hilo  data  of  Cox.    The  Gumbel  distribution  yields  a  l-in-100-year 
elevation  of  28.8  ft  for  data  from  1837  through  1976  and  an  elevation  of 
^2.5  ft  for  data  from  19*+6  through  1976.   This  compares  with  the  eleva- 
tions of  27.3  and  kk.2   ft  calculated  for  the  same  time  periods  using 
the  logarithmic  distribution.   Clearly,  the  arguments  used  earlier  con- 
cerning the  period  of  time  that  must  be  considered  for  a  valid  analysis 
are  not  dependent  upon  the  assumed  frequency  distribution.   Since  ample 
precedents  exist  for  using  the  logarithmic  distribution  for  analyzing 
tsunami  frequency  of  occurrence,  this  distribution  was  used  in  this 
report . 
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PART  III:   NUMERICAL  MODEL 

Description 

kO.      The  interaction  of  tsunamis  with  the  Hawaiian  Islands  was 
determined  by  using  a  hybrid  finite  element  numerical  model  developed 

recently  for  harbor  oscillation  and  wave  scattering  problems  by  Chen 

23 

and  Mei   at  Massachusetts  Institute  of  Technology.   The  model  solves 

the  following  generalized  HeLmholtz  equation: 

2 
V[d(x,y)V<j>(x,y)]  +  3-  4>U,y)  =  0  (2) 


where 

V  =  gradient  operator,  ft 

x,y  =  Cartesian  coordinates,  ft,  of  the  location 

d(x,y)  =  water  depth  at  the  location 

<f>(x,y)  =  total  velocity  potential  at  the  location  with  U(x,y), 
a  two-dimensional  vector,  equal  to  -V<j)(x,y) 

a)  =  angular  frequency 

g  =  acceleration  due  to  gravity 

Equation  2  governs  small  amplitude  undamped  long  waves  in  a  region  with 
land  masses  of  arbitrary  shape  and  water  of  variable  depth.   It  has 
further  been  assumed  that  the  flow  is  irrotational. 
kl.      The  Helmholtz  equation  expressed  as 

?  n2 

V  <j>(x,y)  +  ^  4>U,y)  =  0  (3) 

is  the  governing  equation  for  a  constant-depth  ocean  region  outside  the 
region  containing  the  islands. 

k2.      Waves  incident  upon  islands  in  an  infinite  ocean  produce  a 
scattered  wave  having  a  velocity  potential  <\>        given  by 

o 
oo 

4>0  =  >  H  (kr)(a  cos  n  9  +  3  sin  n  6)  {h) 

o         Z-<   n      n  n 

n=0 
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where 

H  =  Hankel  function  of  the  first  kind  of  order  n 
n 

k  =  wave  number 

r  =  spherical  coordinate,  ft 

a  =  unknown  coefficient 
n 

6  =  spherical  coordinate,  radians 

3  =  unknown  coefficient 
n 

k3.      The  symbol   <j>   satisfies  the  Sommerfeld  radiation  condition 

o 

that  the  scattered  wave  must  behave  as  an  outgoing  wave  at  infinity  and 
may  be  expressed  mathematically  as 


lim  /T(  |-  -  ik  )  <f>c  =  0 


(5) 


where   i   equals   v-1  . 

kk.      Chen  and  Mei  used  a  calculus  of  variations  approach  and  ob- 
tained a  Euler-Lagrange  formulation  of  the  boundary  value  problem.   The 
following  functional  with  the  property  that  it  is  stationary  with  re- 
pect  to  arbitrary  first  variations  of  <j>(x,y)  was  constructed  by  Chen 
and  Mei : 
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where 


R  =  region  containing  the  islands 

§  =  line  integral 

$_,  =  far-field  velocity  potential 

<j)  =  velocity  potential  of  incident  wave 

n  =  unit  normal  vector  outward  from  region  R 
a  D 

a  =  boundary  of  region  R 
(J)  =  total  velocity  potential  evaluated  on  boundary  a 

Si 
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1+5.   Proof  was  given  by  Chen  and  Mei  that  the  stationarity  of  this 
functional  is  equivalent  to  the  original  boundary  value  problem. 

k6.      The  integral  equation  obtained  from  extremizing  the  functional 
is  solved  by  using  the  finite  element  method.   This  method  is  a  technique 
of  numerical  approximation  that  involves  dividing  a  domain  into  a  number 
of  nonoverlapping  subdomains ,  which  are  called  elements. 

hj .      The  solution  of  the  problem  is  approximated  within  each  ele- 
ment by  suitable  interpolation  functions  in  terms  of  a  finite  number  of 
unknown  parameters.   These  unknown  parameters  are  the  values  of  the 
field  variable   c|)(x,y)   at  a  finite  number  of  points,  which  are  called 
nodes.   The  relations  for  individual  elements  are  combined  into  a  system 
of  equations  for  all  unknown  parameters. 

U8.   In  the  constant  depth  region  outside  the  region  containing 
the  islands,  the  velocity  potentials  are  solved  analytically  in  terms 
of  unknown  coefficients.   The  region  is  considered  to  be  a  single  ele- 
ment with  an  "interpolation  function"  given  by  Equation  h.      The  infinite 
series  is  terminated  at  some  finite  value  such  that  the  addition  of 
further  terms  does  not  significantly  influence  the  calculated  values  of 
<j)(x,y)  .   The  resulting  equation  is  combined  with  the  system  of  equa- 
tions for  unknown  parameters  at  nodal  points  within  the  region  contain- 
ing the  islands,  and  this  complete  system  is  solved  using  Gaussian 
elimination  matrix  methods. 

1+9.   The  relation  of  n(x,y)   to   <f>(x,y)   through  the  linearized 
dynamic  free  surface  boundary  condition  is  expressed  as 


n(x,y)  -  -  — 


94>(x,y) 

at 


(7) 


where 

n  =  wave  amplitude 

t  =  time 

50.   The  hybrid  finite  element  method  (so  named  by  Chen  and  Mei 
because  the  method  involves  the  combination  of  analytical  and  finite 
element  numerical  solutions)  is  a  steady-state  solution  of  the  boundary 
value  problem.   The  response  of  a  group  of  islands  to  an  arbitrary 
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tsunami  can  be  easily  determined  within  the  framework  of  a  linearized 
theory.   For  example,  an  arbitrary  tsunami  in  the  deep  ocean  can  be 
Fourier  decomposed  as  follows : 


oo 

/ 


b0(t  )e-i[«t+p(»)]  , 


where 


b  =  incident  wave  amplitude 
o 


b(to)  =  amplitude  of  frequency  component   to 

p  ( uj )  =  phase  angle 
If   n(x,y,to)   is  the  response  amplitude  at  any  point  (x,y)  along  the 
island  coasts  due  to  an  incident  plane  wave  of  unit  amplitude  and  fre- 
quency co  ,  then  the  response  of  the  islands  to  the  arbitrary  tsunami 


time  history  b  (t)   is  given  by 


C(x,y,t)  =  r)  rb(oJ)n(x,y,03)e-i[wt+p(w)]  du/ 


(9) 


where  the  operation  R  {}  takes  the  real  part  of  the  braced  quantity. 
Therefore,  as  soon  as   n(x,y,to)   is  known  for  all  to  ,  the  island 
response  to  an  arbitrary  tsunami  can  be  calculated. 

51.  Chen  and  Mei  verified  their  finite  element  numerical  model 
by  comparing  the  model's  calculations  with  analytical  and  experimental 
results  for  simple  wave  problems.   Figure  2  shows  a  comparison  of  ana- 
lytical and  numerous  results  of  wave  diffraction  off  a  vertical  cir- 
cular cylinder. 

52.  The  computer  program  of  Reference  23  "was  modified  before 
calculations  were  made  for  tsunami  response  of  the  Hawaiian  Islands. 
Slight  modifications  were  necessary  to  allow  the  model  to  handle  problems 
with  variable  depths,  since  Chen  and  Mei  had  applied  their  original 
numerical  model  only  to  constant  depth  problems.   The  subroutine  that 
solves  the  large  system  of  algebraic  equations  arising  from  finite 
elements  was  modified  to  be  more  efficient  by  taking  advantage  of  the 
sparseness  (many  zero  terms)  of  the  system  of  equations.   This  modifi- 
cation resulted  in  a  decrease  in  the  computational  time  of  the  numerical 
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Figure  2.   Comparison  of  numerical 
and  analytical  results  for  wave 
scattering  off  a  vertical  circu- 
lar cylinder^ 
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model  by  a  factor  of  approximately  six  for  the  grid  shown  in  the  next 
section.   The  solution  of  the  system  of  algebraic  equations  also  had  to 
be  modified  so  that  calculations  involved  only  small  blocks  of  equations 
at  a  time  since  there  were  more  than  800,000  terms  in  the  system  of 
equations. 

Numerical  Grid 

53.   The  finite  element  model  has  great  advantages  relative  to 
finite  difference  numerical  models  in  accuracy  representing  land  shapes 
and  ocean  bathymetry  since  finite  element  techniques  allow  dramatic 
changes  in  element  sizes  and  shapes.   Element  sizes  can  be  large  in  deep 
water  because  wavelengths  are  large  and,  therefore,  adequately  resolved 
by  a  crude  grid.   As  a  wave  enters  shallower  water,  its  wavelength  de- 
creases and  the  elements  of  the  grid  can  be  telescoped  to  smaller  sizes 
with  no  loss  of  resolution. 
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5U.   The  finite  element  grid  in  Figure  3  shows  elements  of  the  grid 
telescoping  from  large  sizes  in  the  deep  ocean  to  triangles  with  areas 

as  small  as  0.5  sq  km  in  shallow  coastal  waters.   For  comparison,  a 

2k 

recent  study  by  Bernard   used  a  finite  difference  grid  that  included 

the  major  islands  of  Hawaii  and  had  grid  cells  with  areas  of  approxi- 
mately 30  sq  km  (square  cells  5.5  km  on  a  side).   Since  the  islands  of 
Hawaii  have  small  surface  areas  and  are  surrounded  by  waters  of  rapidly 


Figure  3.   Finite  element  grid 
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varying  bathymetry,  an  accurate  representation  of  the  islands  in  a  nu- 
merical model  requires  the  small  grid  units  that  are  feasible  (from  an 
economic  viewpoint)  for  a  finite  element  model  only. 

55.  The  geometric  shapes  of  the  eight  islands  comprising  the 
Hawaiian  Island  chain  are  obviously  modeled  very  precisely  by  the  finite 
element  grid  (Figure  3).   Extremely  rapid  bathymetry  changes  also  are 
very  accurately  modeled  by  the  grid.   Furthermore,  the  number  of  node 
points  (corner  points  of  the  triangles)  along  the  shoreline  of  the 
islands  is  very  dense.   Since  wave  amplitudes  are  calculated  at  node 
points  by  the  finite  element  model,  the  model  can  adequately  resolve  the 
rapid  wave  height  variations  along  coastlines  that  are  known  to  occur 
during  tsunami  activity. 

56.  The  finite  element  model  calculations  for  the  grid  (Figure  3) 
are  extremely  rapid.   The  main  calculation  time  involves  the  solution  of 
a  very  large  matrix  with  over  800,000  terms,  yet  the  computational  time 
for  the  model  to  solve  the  response  of  the  islands  to  waves  of  a  partic- 
ular period  is  less  than  1  min  of  computer  time  on  a  Control  Data  Corpo- 
ration 76OO.   This  rapid  calculation  time  makes  it  economically  feasible 
to  determine  the  response  of  the  islands  to  arbitrary  tsunamis,  as  dis- 
cussed in  the  following  section.   One  reason  that  the  computational 
time  required  by  the  finite  element  model  is  so  much  less  than  that  re- 
quired by  finite  difference  models  is  that  the  finite  element  model  uses 
small  elements  only  in  areas  where  they  are  necessary.   The  grid  in 
Figure  3  has  approximately  2,500  nodal  points,  whereas  the  finite  dif- 
ference  grid  of  Bernard   had  26,000  grid  points.   Even  so,  some  of  the 
elements  of  the  finite  element  grid  are  up  to  60  times  smaller  than  the 
finite  difference  grid  cells. 

Verification 

57-   The  finite  element  numerical  model  is  verified  in  this  study 
by  comparing  numerical  model  calculations  with  tide  gage  recordings  of 
the  i960  and  196h   tsunamis.   The  i960  Chilean  and  196k   Alaskan  tsunamis 
are  the  only  major  tsunamis  for  which  some  reliable  information 


31 


regarding  source-generating  characteristics  exists,  with  much  more  in- 
formation existing  for  the  Alaskan  source  than  the  Chilean  source. 

58.   A  deepwater  recording  of  a  tsunami  far  from  the  perturbing 
influences  of  a  coastal  area  has  never  been  made.   Therefore,  a  proto- 
type wave  record  in  deep  water  of  the  i960  Chilean  and  196I*  Alaskan 
tsunamis  is  not  available  for  use  as  input  to  verify  the  finite  element 

model.   However,  finite  difference  numerical  models  have  been  used  in 

11  25  26  27 
previous  studies   'to  simulate  the  uplift  deformation  of  the 

ocean  water  surface  caused  by  the  permanent  vertical  displacement  of 
the  ocean  bottom  during  an  earthquake  and  the  subsequent  propagation  of 
the  tsunami  across  the  deep  ocean.   The  permanent  deformation  (permanent 
in  the  sense  that  the  time  scale  associated  with  it  is  much  longer  than 
the  period  of  the  tsunami)  is  considered  to  be  the  important  parameter 
in  tsunami  generation  and  not  the  transient  movements  within  the  time 
history  of  the  ground  motion.   These  transient  movements  occur  over 
periods  of  time  of  the  order  of  seconds,  whereas  tsunami  wave  periods 
are  of  the  order  of  tens  of  minutes.   Reference  28  confirms  the  fact 
that  the  permanent  ground  displacement  of  the  source  region  and  not  the 
transient  ground  movements  determines  the  characteristics  of  the  result- 
ing tsunami. 

59-   Figures  h   and  5  show  time  histories  of  the  i960  Chilean  and 

196h   Alaskan  tsunamis  in  deep  water  near  the  Hawaiian  Islands  calculated 

11  27 
by  finite  difference  models  in  previous  studies.   '    The  permanent 

deformation  of  the  ocean's  bottom  at  the  source  as  a  function  of  spatial 
location  was  taken  from  Reference  29  for  the  Alaskan  source  and  Refer- 
ences 30  and  31  for  the  Chilean  source.   Grids  measuring  1/3  by  1/3  deg 
were  used  for  the  calculations. 

60.   The  hybrid  finite  element  method  is  a  steady-state  solution 
of  the  boundary  value  problem.   As  described  in  paragraphs  1*0-52,  the 
response  of  a  location  in  the  Hawaiian  Islands  to  any  tsunami  can  be 
calculated  by  the  finite  element  model  once  the  response  amplitude 
n(x,y,a))   due  to  an  incident  plane  wave  of  unit  amplitude  and  frequency 
03  is  known  at  the  location   (x,y)   for  all  frequencies.   Of  course,  it 
is  not  feasible  to  calculate  the  integrals  of  Equations  8  and  9  over  all 
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frequencies.   Instead  the  frequency  range  is  discretized,  and  the  in- 
tegrals replaced  by  sums  over  a  frequency  range  containing  most  of  the 
energy  of  the  tsunami. 

61.  Equation  8  involves  a  Fourier  decomposition  of  a  time  series. 
This  decomposition  vas  accomplished  for  the  time  histories  of  the  i960 
Chilean  and  I96U  Alaskan  tsunamis  (Figures  k   and  5)  by  using  a  least- 
squares  harmonic  fitting  procedure.   The  time  history  of  the  Alaskan 
tsunami  in  deep  water  was  decomposed  into  18  components  with  periods 
ranging  from  lU.5  min  to  the  time  length  of  the  record.   The  variance 
of  the  residual  (difference  between  the  actual  record  and  a  recomposi- 
tion  of  the  18  components)  was  approximately  0.2  percent  of  the  variance 
of  the  record.   Therefore,  virtually  all  the  energy  in  the  wave  record 
was  contained  in  the  18  components.   The  time  history  of  the  Chilean 
tsunami  in  deep  water  was  decomposed  into  11  components  with  periods 
ranging  from  15-5  min  to  the  time  length  of  the  record.   The  variance 

of  the  residual  was  less  than  0.1  percent  of  the  variance  of  the  record. 
For  both  cases,  the  original  time  history  and  a  time  history  constructed 
from  a  recomposition  of  the  components  were  virtually  indistinguishable 
visually. 

62.  Equations  8  and  9  take  the  following  form  when  discretized: 

A       -i[w  t+p(u>  )] 
b  (t)  =  >  b(u  )e    n     n  (10) 

n=i 


and 


\  A  -i[o)  t+p(u  )]  J 

C(x,y,t)  =  Re  O  b(u>n)n(x,y,u)n)e    n     n  \  (ll) 

where 

m  =  number  of  components 

go  =  frequency  of   nth  components 

63.   The  b(u)  )   term  is  determined  by  the  least-squares  harmonic 
fitting  procedure,  and  n(x,y,oj  )   is  determined  by  the  finite  element 
numerical  model  for  each  frequency  oj   and  at  each  location  (x,y). 
Therefore,  the  time  history   £(x,y,t)  ,  which  represents  the  response 
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of  location  (x,y)  to  the  deepwater  tsunami  time  history  represented  by 
b  (t)  ,  can  be  simply  calculated  for  any  location  along  the  coastline 
of  the  Hawaiian  Islands. 

6k.      Figures  6  through  8  present  a  comparison  between  tide  gage 
recordings  of  the  196k   tsunami  at  Kahului  (Maui),  Honolulu  (Oahu),  and 
Hilo  (Hawaii),  respectively,  and  the  numerical  model  calculations.   The 
largest  waves  recorded  at  each  of  these  sites  are  shown,  but  the  smaller 
waves  arriving  at  later  times  and  thus  of  little  interest  to  this  study 
are  not.   Whenever  gage  limits  were  encountered,  the  tide  gage  record- 
ings were  linearly  extended.   Of  course,  a  tide  gage  is  a  nonlinear 
device  with  a  response  depending  both  upon  the  period  and  amplitude  of 
the  disturbance  measured.   However,  based  upon  the  work  of  Noye ,   the 
distortion  of  the  tsunami  shown  in  Figures  6  through  9  by  a  standard 
tide  gage  is  not  significant. 

65.  The  comparisons  shown  in  Figures  6  through  8  are  in  remarkable 
agreement,  especially  considering  the  fact  that  the  ground  displacement 
of  the  196U  earthquake  was  not  precisely  known.   Since  the  Hilo  break- 
water was  not  included  in  the  numerical  model  grid,  the  numerical  model 
calculations  for  Hilo  are  too  large.   This  breakwater  was  undoubtedly 
highly  permeable  to  the  196k   tsunami.   However,  it  probably  did  cause 
some  energy  dissipation  and  reflection,  although  how  much  is  not  known. 
The  numerical  model  results  appear  to  be  too  large  in  Hilo  by  some  con- 
stant factor,  since  the  tide  gage  recording  and  the  numerical  model  cal- 
culations have  the  same  form  with  the  first  wave  crest  and  trough  having 
approximately  the  same  amplitude  and  being  proportionately  smaller  than 
the  second  crest. 

66.  Figure  9  shows  a  comparison  between  the  tide  gage  recording 
of  the  i960  tsunami  at  Honolulu  and  the  numerical  model  calculations. 
The  Honolulu  gage  was  the  only  tide  gage  in  the  Hawaiian  Islands  not 
seriously  damaged  by  this  tsunami.   Again  the  largest  waves  recorded 
are  shown.   The  agreement  indicated  in  Figure  9  is  quite  good  consider- 
ing the  fact  that  the  ground  displacement  of  the  i960  earthquake  was 
not  known  precisely.   Figure  9  is  the  first  comparison  ever  published 
of  a  recording  of  the  i960  tsunami  and  a  numerical  model  calculation. 
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Figure  6.   196^  tsunami  from  Alaska  recorded  at  Kahului ,  Maui 
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Figure  7.   196^  tsunami  from  Alaska  recorded  at  Honolulu,  Oahu 
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Figure  8.   1964  tsunami  from  Alaska  recorded  at  Hilo,  Hawaii 
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67.  The  comparisons  displayed  in  Figures  6  through  9  are  also  the 
first  close  comparisons  of  tide  gage  recordings  of  tsunamis  in  shallow 
water  and  numerical  model  calculations.   Numerical  model  calculations 
have  been  compared  with  a  recording  of  the  196U  tsunami  by  Van  Dora's 

gage  at  Wake  Island  which  was  located  at  a  depth  of  approximately 

33 
800  ft.     Only  the  first  wave  of  the  tsunami  was  in  good  agreement  in 

the  comparison,  however.   A  comparison  also  has  been  previously  made  for 
the  I96U  tsunami  at  Hilo.     However,  this  comparison  took  a  time  history 
of  the  196h   tsunami  calculated  over  100  miles  from  Hilo  in  a  water  depth 
of  approximately  15,000  ft  and  used  it  as  input  to  a  fine  finite  dif- 
ference grid  covering  only  the  immediate  Hilo  area  with  an  input  bound- 
ary approximately  6  to  7  miles  from  Hilo  in  about  1,000  ft  of  water. 
The  wave  period  and  phase  of  the  waves  in  the  comparison  were  not  in 
good  agreement ,  but  the  amplitudes  of  the  largest  waves  were  in  reason- 
able agreement.   The  amplitudes  of  the  computed  waves  depend  upon  the 
water  depth  at  the  input  boundary  of  the  grid,  however,  and  can  be 
changed  by  varying  this  depth.   No  explanation  was  given  in  Reference  3^ 
for  the  depth  used;  however,  it  seems  likely  that  it  was  chosen  because 
the  computed  waves  had  reasonable  amplitudes  when  this  depth  at  the  in- 
put boundary  was  tried  and  not  when  others  were  tried. 

Model  Use 

68.  The  use  of  a  numerical  model  as  an  aid  for  interpolation 
between  historical  data  points  can  be  illustrated  for  the  i960  tsunami 
from  Chile  (there  were  too  few  measurements  of  the  I96U  tsunami  from 
Alaska  to  be  illustrative).   Figure  10  shows  a  comparison  on  Oahu  of  the 
shoreline  elevations  calculated  by  the  numerical  model  with  historical 
measurements  for  the  i960  tsunami.   The  numerical  and  historical  eleva- 
tions do  not  compare  as  favorably  as  the  comparison  shown  in  Figure  9 
for  several  reasons.   First,  the  historical  measurements  are  not  at  the 
shoreline.   As  discussed  in  Part  II,  these  measurements  occurred  on  an 
average  of  200  ft  from  the  shoreline.   The  water  depths  in  elements  of 
the  numerical  grid  touching  land  also  are  variable  around  the  island. 
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The  depths  cannot  approach  zero  and,  therefore,  the  shoaling  process  may 
not  he  complete.   The  tide  gage  recording  shown  in  Figure  9  was  in  a 
depth  of  water  well  represented  in  the  numerical  grid.   Furthermore, 
some  parts  of  Oahu  have  reefs  that  may  significantly  influence  tsunamis 
and  may  not  be  well  simulated  by  the  numerical  model  (especially  if  they 
are  relatively  narrow).   Only  the  effect  on  the  i960  tsunami  of  the 

reefs  of  Kaneohe  Bay  are  known  from  historical  data.   A  tide  gage  at 

35 
Mokuoloe  Island  in  Kaneohe  Bay   clearly  shows  that  the  reefs  of  the  bay 

act  like  a  low-pass  filter,  filtering  out  wave  periods  below  approxi- 
mately 1  hr.   Thus,  only  components  with  periods  above  1  hr  were  recom- 
posed  to  give  the  elevations  in  Kaneohe  Bay  shown  in  Figures  10 
through  16. 

69.  Figure  11  shows  a  comparison  of  the  numerical  model  eleva- 
tions of  Figure  10  multiplied  by  correction  factors  with  historical 
measurements  for  the  i960  tsunami  from  Chile.   The  correction  factor  for 
a  point  is  determined  by  taking  the  ratio  of  the  nearest  historical  mea- 
surement elevation  and  the  numerical  model  elevation  at  the  historical 
measurement  location.   This  correction  factor  is  then  multiplied  by  the 
numerical  model  elevation  at  the  point  of  interest.   The  final  step  is  to 
require  that  the  curve  through  the  corrected  numerical  model  elevations 
also  pass  through  any  historical  measurements  through  which  it  does  not 
already  pass.   Figure  12  shows  elevations  of  the  i960  tsunami  from  Chile 
determined  by  this  process.   The  correction  factor  determines  all  of  the 
local  factors  not  already  determined  by  the  numerical  model.   These 
local  factors  include  the  effects  of  land  flooding  to  the  point  of  the 
historical  measurement,  final  shoaling  near  the  shoreline,  and  effects 

of  reefs.   Figure  12  is,  therefore,  a  reconstruction  of  elevations  for 
the  i960  tsunami  from  Chile  on  the  island  of  Oahu. 

70.  Figure  13  shows  a  comparison  on  Oahu  of  normalized  numerical 
model  shoreline  elevations  for  a  generalized  tsunami  from  Kamchatka  with 
normalized  historical  measurements  for  the  1952  tsunami  from  Kamchatka. 
Since  the  vertical  uplift  of  the  ground  in  the  source  region  is  not 
known  for  this  tsunami,  the  deepwater  form  of  the  tsunami  cannot  be  de- 
termined using  a  numerical  model,  as  was  done  for  the  i960  and  I96U 
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tsunamis.   Therefore,  an  average  or  generalized  Kamchatkan  tsunami  must 
be  used  as  input  to  the  numerical  model,  as  discussed  in  Part  II.   Since 
a  deep  ocean  height  of  the  1952  tsunami  or  a  generalized  Kamchatkan 
tsunami  is  not  known,  normalized  elevations  are  shown.   Both  the  numeri- 
cal model  elevations  and  the  historical  data  have  been  normalized  on 
each  of  the  four  sides  of  Oahu  by  the  largest  corresponding  elevation  on 
that  particular  side.   Figures  ik   and  15  show  comparisons  for  the  Kam- 
chatkan tsunami  analogous  to  Figures  11  and  12  for  the  i960  tsunami. 
The  normalization  has  been  eliminated  in  Figure  15  by  multiplying  the 
elevations  on  each  island  side  by  the  inverse  of  the  factor  used  to 
normalize  the  historical  data.   Figure  16  shows  elevations  of  the  I9U6 
tsunami  from  the  Aleutians  determined  from  the  historical  data  and  the 
numerical  model  calculations.   Again,  a  generalized  tsunami  had  to  be 
used  to  represent  the  19^+6  tsunami  because  the  source  characteristics 
and  consequently  the  deepwater  form  of  this  tsunami  are  unknown. 

71.  Figures  12,  15,  and  l6  can  be  used  to  determine  elevations  on 
Oahu  for  the  i960,  1952,  and  19^6  tsunamis,  respectively,  even  for  loca- 
tions not  having  historical  data  for  any  or  all  of  these  tsunamis.  Sim- 
ilar curves  for  the  1957  and  196^+  tsunamis  were  derived.  The  same  tech- 
niques were  used  to  allow  the  prediction  of  elevations  for  these  five 
tsunamis  on  all  of  the  islands  of  Hawaii. 

72.  Reference  36  also  used  a  numerical  model  (the  finite  differ- 
ence  model  developed  by  Bernard   )  to  interpolate  between  estimated 
(from  historical  data)  elevations  produced  on  five  islands  by  a  tsunami 
from  the  north  with  a  magnitude  of  four.   This  finite  difference  model 
had  a  grid  size  much  too  crude  to  adequately  resolve  island  shapes, 
ocean  bathymetry,  or  tsunami  wavelengths  (see  paragraphs  5^-56).   A 
single  Gaussian  wave  crest  is  inputted  into  the  model.   If  a  resolution 
of  four  grid  cells  per  crest  wavelength  is  maintained,  this  model  cannot 
propagate  waves  in  water  with  a  depth  less  than  about  500  ft.   The  main 
processes  involved  in  the  transformation  of  the  form  of  a  tsunami  prob- 
ably occur  in  the  region  over  the  continental  shelf  (depths  less  than 
600  ft).   Consequently,  this  model  calculates  elevations  much  smaller 
than  historical  data.   The  effects  of  propagation  over  the  continental 
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shelf  are  determined  by  multiplying  the  numerical  model  elevations  by  a 
single  factor  for  each  island  (two  factors  for  the  island  of  Hawaii), 
which  results  in  the  best  agreement  between  the  historical  data  and  the 
corrected  numerical  model  elevations. 

73.   The  agreement  is  poor  between  historical  data  and  the  cor- 
rected numerical  model  elevations  presented  in  Reference  36.   For  exam- 
ple, on  the  northeast  side  of  Oahu,  from  Kahuku  Point  to  Makapuu  Point, 
the  elevations  for  a  magnitude  four  tsunami  from  the  north  vary  from 
7  to  10.7  ft.   This  remarkable  lack  of  variation  is  undoubtedly  a  result 
of  using  a  single  factor  for  Oahu  to  represent  the  effects  of  propaga- 
tion over  the  continental  shelf.   For  a  magnitude  five  tsunami,  such  as 
the  19^6  Aleutian  tsunami,  these  elevations  must  be  multiplied  by  a 
factor  of  two,  resulting  in  an  elevation  range  from  ik   to  21. k   ft. 
Reference  10  presents  36  historical  elevations  for  the  19^+6  Aleutian 
tsunami  on  the  northeast  coast  of  Oahu  from  Kahuku  Point  to  Makapuu 
Point,  and  only  7  fall  within  this  range  of  elevations.   The  historical 
elevations  actually  vary  from  1  to  37  ft  on  this  coast. 
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PART  IV:   RESULTS 


Use  of  Plots 


jh.      Plates  1-8  show  the  location  numbers  that  are  plotted  along 
the  location  axis  in  Plates  9-^-   Plate  2  shows  the  location  numbers 
for  the  Hilo  area.   The  location  numbers  are  equidistantly  plotted  in 
Plates  9-^  even  though  the  distance  between  them  shown  in  Plates  1-8 
is  not  equidistant.  Plates  9-443    therefore  show  location  and  not  dis- 
tance along  the  coast.      However^   distances  are  true  between  integer  num- 
bers plotted  along  the   location  axis.      That  is,  if  a  place  of  interest 
is  located  halfway  between  two  location  numbers  shown  on  one  of  the 
plots  of  Plates  1-8,  it  is  also  located  halfway  between  the  correspond- 
ing location  numbers  of  Plates  9-kh. 

75.   Plots  displaying  the  A  and  B  coefficients  are  shown  in 
Plates  9-^«   An  island  by  island  (plus  the  Hilo  area)  breakdown  of  the 
A-  and  B-coefficient  plots  is  as  follows  (the  uninhabited  U.  S.  Navy 
target  island  of  Kahoolawe  (Figure  l)  is  not  included): 


Island 

Coefficient  A 

Coefficient  B 

Hawaii 

Plates  9-lU 

Plates  15-20 

(Hilo) 

Plates  12 

Plate  18 

Oahu 

Plates  21-23 

Plates  2U-26 

Maui 

Plates  27-29 

Plates  30-32 

Kauai 

Plates  33-31* 

Plates  35-36 

Molokai 

Plates  37-38 

Plates  39-^0 

Lanai 

Plate  kl 

Plate  \2 

Niihau 

Plate  U3 

Plate  \h 

76.   To  illustrate  how  A  and  B  are  determined  using 
Plates  1-kk ,   the  city  of  Waiahole  in  Kaneohe  Bay  on  the  island  of  Oahu 
has  a  location  number  of  31  (Plate  3).   Plate  21  gives  a  value  of  2.0  for 
coefficient  A  ,  and  Plate  2k   yields  a  B  coefficient  of  1.5-   Then, 
these  values  of  A  and  B  plus  F  ,  which  equals  0.05  for  a  l-in-20- 
year  event,  are  substituted  into  Equation  1  (see  paragraph  37)  as  follows: 


h20  =  ~1,5  "  2'°   logio(0-05) 
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=  -1.5  -  2.0(-1.3) 
=  1.1  ft 

77 •   If  F  equals  0.01  for  a  l-in-100-year  event,  substituting 
in  Equation  1  results  in 


=  -1.5  -  2.0(-2) 
=  2.5  ft 


78.   The  tsunami  elevations  in  Kaneohe  Bay  are  very  small  as  a 
result  of  an  extensive  reef  protecting  the  hay  from  tsunamis. 

79-   Another  example  is  a  calculation  for  Makawana  Point  on  the 
island  of  Maui.   Plate  h   shows  Makawana  Point  lying  one  third  of  the 
way  between  points  3^  and  35 •   The  values,   A  =  22  and  B  =  22.5  >  come 
from  Plates  28  and  31,  respectively.   Substituting  these  values  in 
Equation  1  with  a  l-in-20-year  frequency  yields 


h2Q  =  -22.5  -  22  log  0.05 
=  -22.5  -  (22)(-1.3) 
6.1  ft 


and  with  a  l-in-100-year  frequency  gives 


h10Q  =  -22.5  -  22  log  0.01 
=  -22.5  -  (22)(-2) 
=   21.5  ft 


One-in-Ten-Year  Heights 

80.   Part  II  discusses  the  fact  that  the  linear  relation  between 
the  logarithm  of  tsunami  frequency  of  occurrence  and  tsunami  elevations 
at  Hilo  is  valid  for  frequencies  of  occurrence  lower  than  0.1  per  year 
only.   There  are  some  locations  where  this  linear  relation  is  valid  for 
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frequencies  of  occurrence  equal  to  or  lower  than  0.05  per  year  only.   It 
is,  therefore,  recommended  that  Equation  1  and  Plates  9-U4  be  used  for 
frequencies  of  occurrence  equal  to  or  lower  than  0.05  per  year.   The  one- 
in-ten-year  frequency  of  occurrence  for  the  entire  Hawaiian  coastline 
can  be  determined  by  using  Table  1.   The  elevations  presented  in  this 
table  were  estimated  by  multiplying  the  elevation  of  the  tenth  largest 
tsunami  at  a  location  (return  period  of  lUO  years  divided  by  ten,  i.e. 
lU  years)  by  0.7-   The  value  of  0.7  is  the  ratio  of  the  return  period  of 
the  frequency  0.1  tsunami  (10  years)  and  the  return  period  of  the  tenth 
largest  tsunami  (l^+  years).   Frequencies  of  occurrence  between  0.1  and 
0.05  can  be  linearly  interpolated  between  the  0.1  values  from  Table  1 
and  calculations  using  Equation  1  for  a  l-in-20-year  frequency  of 
occurrence. 

Discussion 

81.  Since  the  plots  in  Plates  9-^  are  strongly  based  upon  histor- 
ical data,  they  should  show  trends  similar  to  those  in  the  historical 
data.   Of  course,  frequency  of  occurrence  depends  upon  both  the  absolute 
and  relative  magnitudes  of  the  A  and  B  coefficients.   It  is  difficult 
to  make  generalizations  involving  only  one  of  these  coefficients.   How- 
ever, it  is  true  that  the  elevations  calculated  for  tsunamis  with  a  low 
frequency  of  occurrence  are  usually  large  when  A  is  large. 

82.  Plates  9-1^-  show  the  A  coefficient  for  the  big  island  of 
Hawaii.   The  coefficient  is  fairly  small  all  along  the  west  coast  of 
Hawaii  except  for  moderate  values  from  Napoopoo  to  Keauhau  resulting 
from  the  fairly  large  tsunami  that  struck  this  area  in  1896.   However, 
A  is  large  along  the  northeast  coast  of  Hawaii  with  a  very  large  peak 
near  the  Pololu  Valley  caused  by  the  large  tsunami  elevations  that 
occurred  during  the  19^6  and  i960  tsunamis.   In  the  Hilo  Bay  area,  the 

A  value  is  fairly  large  but  not  in  comparison  with  the  A  values  along 
the  northeast  coast.   Extensive  damage  and  great  loss  of  life  occurs  in 
Hilo  but  not  along  the  northeast  coast  because  Hilo  is  the  most  populous 
city  on  the  island  of  Hawaii  and  has  a  low- lying  waterfront  area, 
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whereas  the  northeast  coast  is  sparsely  settled  with  villages  on  high 
cliffs.   Plates  13  and  lU  show  large  A  values  on  the  South  Puna 
coast  due  to  the  local  tsunamis  of  1868  and  1975 • 

83.   Plates  21-23  show  A  values  for  the  island  of  Oahu.   The 
largest  A  values  for  Oahu  are  generally  at  the  corners  or  tips  of  the 
island.   For  example,  there  are  peaks  in  the  A-coefficient  plots  at 
Kaena  Point,  Kahuku  Point,  Mahie  Point,  Mokapu  Point,  Makupuu  Point, 
and  Diamond  Head  (Figure  l).   Other  islands  of  the  Hawaiian  Islands 
often  display  a  similar  trend  of  large  A  values  on  island  tips.   The 
A  coefficient  values  are  small  on  Oahu  in  reef-protected  Kaneohe  Bay 
and  in  Pearl  Harbor. 

Qk.      The  l-in-100-year  elevations  for  Oahu  that  can  be  calculated 
using  the  information  presented  in  this  report  are  in  general  less  than 
the  l-in-100-year  elevations  calculated  in  an  earlier  tsunami  study 
(Reference  k)    except  at  Honolulu.   However,  Reference  k   used  data  for 
the  period  19^6  through  1964  only,  except  for  Honolulu,  where  data 
dating  back  to  1837  were  used. 

85.  Values  of  A  for  the  island  of  Maui  are  shown  in  Plates  27- 
29.   The  A  coefficient  is  fairly  small  on  the  south  coast  of  the 
island  and  large  on  the  north  coast.   This  difference  in  A  values 
between  the  north  and  south  coasts  of  the  island  is  due  to  the  fact  that 
the  I9J+6  tsunami  from  the  Aleutian  Islands  produced  the  largest  histori- 
cal elevations  on  the  island.   The  north  coast  is  directly  exposed  to 
tsunamis  generated  in  the  Aleutians,  whereas  the  south  coast  is  shielded. 
This  same  effect  of  wave  elevations  on  the  back  side  of  an  island  being 
less  than  elevations  on  the  side  of  the  island  facing  the  tsunami- 
generating  area  occurs  on  Kauai,  Molokai ,  and  Niihau.   The  19^-6  and 
1957  tsunamis  from  the  Aleutian  Islands  were  the  largest  to  strike  all 
three  of  these  islands.   The  south  coast  of  Molokai  is  also  reef  pro- 
tected, and  this  partially  accounts  for  its  exceptionally  low  values  of 
A. 

86.  Since  data  were  taken  from  a  lUO-year  period  of  time 
(1837-1976),  the  1-in-lUO-year  tsunami  elevation  at  a  location  should 
be  similar  to  the  largest  value  during  the  lUO-year  period  at  the 
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location.   This  is  usually  the  case,  of  course.   For  example,  the 
1-in-l^O-year  elevation  near  Mooheau  Park  in  Hilo  (Plate  2)  is  approxi- 
mately 28  ft  according  to  Equation  1  and  Plates  12  and  18.   The  largest 
historical  elevation  above  mean  sea  level  was  26  ft  during  the  19^-6 
tsunami.   The  1-in-lUO-year  elevation  at  Lahaina,  Maui,  is  10  ft  ac- 
cording to  Equation  1  and  Plates  27  and  30.   The  largest  historical 
elevation  was  9  ft  during  the  i960  tsunami. 

Risk  Calculation 

87.  The  average  frequency  of  occurrence  F  is  a  mean  exceedance 
frequency,  i.e.  an  average  frequency  per  year  of  tsunamis  of  equal  or 
greater  elevation.   It  is  also  possible  to  calculate  the  chance  of  a 
given  elevation  being  exceeded  in  a  certain  period  of  time.   Such  a 
calculation  is  a  risk  calculation. 

88.  Tsunamis  are  usually  caused  by  earthquakes,  and  earthquakes 
are  often  idealized  as  a  generalized  Poisson  process.   '    Many  investi- 
gators have  assumed  that  tsunamis  also  follow  such  a  stochastic  pro- 

19  39 
cess.   '    The  probability  that  a  tsunami  with  an  average  frequency  of 

occurrence  of  F  is  exceeded  in  D  years,  assuming  that  tsunamis  fol- 
low a  Poisson  process,  is  given  by  the  following  equation: 

P  =  1  -  e_FI)  (12) 

89.  For  example,  the  probability  that  the  l-in-100-year  elevation 
of  21.5  ft  at  Makawana  Point,  Maui  (Plate  k) ,   will  occur  in  a  50-year 
period  is 

p-1..-(0.01)(50) 

-  1  -  a"0'05 

=   1   -   0.6l 
=   0.39 

and  in  a  10-year  period  is 
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p  =  i  _  e-(o.oi)(io) 

_  ,    -o.i 

=  1  -  e 
=1-0.9 
=  0.1 

90.  As  mentioned  in  the  introductory  section  of  Part  I,  an  eval- 
uation of  risk  is  important  whenever  human  life  may  "be  exposed  to  pos- 
sible danger  as  a  result  of  land  development.   For  example,  it  would 
certainly  be  extremely  foolhardy  to  base  tsunami  evacuation  zones  simply 
on  something  like  the  l-in-100-year  elevations  given  in  this  report.   As 
shown  in  the  previous  paragraph,  the  odds  that  a  l-in-100-year  event 
will  occur  during  a  10-year  period  are  not  insignificant. 

91.  Risk  calculations  can  be  used  to  add  a  safety  factor  when 
evaluating  whether  or  not  it  is  prudent  to  develop  land  at  some  eleva- 
tion that  might  expose  human  life  to  possible  danger.   First,  it  is 
necessary  to  decide  on  some  acceptable  risk.   For  example,  perhaps  a  one 
chance  in  10,000  that  an  elevation  be  exceeded  during  a  10-year  period 
is  an  acceptable  risk.   (This  risk  is  used  only  for  illustration  and  not 
to  suggest  that  such  a  risk  is  or  is  not  acceptable.)   Then  D  of 
Equation  12  equals  10  and  P  equals  0.0001.   Substituting  these  values 
of  D  and  P  in  Equation  12  yields  an  F  equal  to  approximately 
0.00001.   Therefore,  the  elevation  for  which  there  is  only  one  chance  in 
10,000  that  the  elevation  will  be  exceeded  during  a  10-year  period  is 
the  l-in-100, 000-year  elevation.   For  Makawa  Point,  Maui,  this  elevation 
would  be  87  ft  (substituting  F  =  0.00001  and  A  and  B  values  from 
Plates  28  and  31  into  Equation  l).   Of  course,  Equation  1  may  not  be 
valid  for  the  l-in-100, 000-year  tsunami,  since  the  linear  relationship 
between  h  and  the  logarithm  of  F  may  not  be  valid  for  arbitrarily 
small  values  of  F  .   However,  there  is  not  sufficient  historical  data 
at  this  time  to  determine  any  other  relationship  between  h  and  F  . 
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PART  V:   CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


92.  The  hybrid  finite  element  numerical  model  was  shown  to  be  an 
excellent  numerical  method  for  determining  the  interaction  of  tsunamis 
with  the  Hawaiian  Islands.   The  comparisons  presented  between  the  numer- 
ical model  calculations  and  historical  tide  gage  recordings  for  the  i960 
and  196k   tsunamis  were  the  first  close  comparisons  for  tsunamis  in  shal- 
low water.   This  accurate  modeling  of  tsunamis  by  a  numerical  model  is  a 
significant  breakthrough  in  tsunami  research.   Future  studies  using  this 
numerical  model  may  develop  a  greater  understanding  of  the  processes  that 
lead  to  large  tsunami  elevations  at  one  location  along  a  coast  and  small 
elevations  at  a  nearby  location.   Understanding  a  problem  is  the  first 
step  toward  alleviating  it.   Furthermore,  sound  floodplain  management 
involves  not  only  understanding  possible  hazards  to  an  area  that  may  be 
developed,  but  also  the  effect  the  development  itself  may  have  on  the 
hazard.   The  finite  element  model  can  be  used  to  test  the  effect  on  tsu- 
nami elevations  of  topographic  or  bathymetric  modifications  in  an  area. 

93.  This  report  presents  a  simple  method  to  determine  the  fre- 
quency of  occurrence  of  tsunami  elevations  above  mean  sea  level  at  the 
shoreline  for  frequencies  equal  to  or  lower  than  l-in-10  years  for  the 
entire  coastline  of  the  Hawaiian  islands  (except  for  the  uninhabited 
U.  S.  Navy  target  island  of  Kahoolawe ) .   For  some  of  the  coastline  of 
the  Hawaiian  Islands,  runup  above  mean  sea  level  can  be  shown  to  nearly 
equal  wave  elevation  above  mean  sea  level  at  the  shoreline  both  for 
tsunamis  producing  water-surface  changes  similar  to  rapidly  rising  tides 
and  for  borelike  waves.   Thus,  the  elevations  calculated  using  the 
methods  presented  also  are  often  runup  elevations.   However,  there  are 
places,  such  as  low- lying  or  estuarine  areas,  or  areas  with  large  ele- 
vations at  the  shoreline,  where  the  wave  elevation  decreases  as  the  wave 
moves  across  the  land  so  that  runup  is  less  than  the  elevation  at  the 
shoreline.   Such  an  effect  probably  occurs  whenever  there  is  signifi- 
cant flooding  during  a  tsunami. 


Recommendations 

94.  Bret Schneider  and  Wybro   have  developed  a  one-dimensional 
(l-D)  approach  to  estimate  flooding  that  neglects  time-dependent  effects 
"but  does  include  dissipative  effects.   Since  the  method  is  l-D,  it 
neglects  flow  convergence  and  divergence;  thus,  it  is  valid  only  for 
straight  coastlines  where  land  elevation  contours  are  parallel  to  the 
coast,  and  there  are  no  natural  or  man-made  barriers  (large  rock  out- 
cropping, buildings,  or  elevated  roads)  or  land  roughness  variations 
along  the  coastline  (e.g.  due  to  elimination  of  homogeneous  natural  veg- 
etation or  topography  by  development  and  replacement  by  partial  grass 
plantings  and  land  grading).   Since  the  method  allows  no  time  dependence, 
it  is  not  valid  for  large  low-lying  areas  where  there  would  be  insuffi- 
cient time  for  complete  flooding  or  areas  where  coastline  restrictions 
(e.g.  harbors)  or  sand  dune  barriers  limit  the  rate  of  flooding  and  thus 
make  the  time  dependence  of  the  flooding  important.   This  method  is  not 
completely  adequate  for  bays,  harbors,  land  points,  developed  areas  such 
as  cities,  large  low- lying  areas,  some  sand-dune  protected  areas,  and 
other  areas  where  there  are  topographical,  roughness,  or  coastline  vari- 
ations.  There  are  many  sections  of  coast  in  the  Hawaiian  Islands  where 
assumptions  of  this  method  are  approximately  satisfied.   The  method  is 
suitable  for  application  on  these  largely  unpopulated  open  coasts. 

95.  There  is  not  a  numerical  model  at  this  time  that  includes 
two-dimensional  and  time-dependent  effects  of  land  flooding  during 

a  tsunami.   However,  the  WES  has  developed  an  implicit  finite  difference 

numerical  model  for  simulating  the  tidal  hydrodynamics  of  an  inlet,  bay, 

kl 
or  harbor.     This  model  includes  variable  bottom  roughness,  nonlinear 

advective  terms  in  the  momentum  equations,  treatment  of  regions  that  are 

inundated  during  a  portion  of  the  tidal  cycle,  exposed  and  submerged 

barriers,  and  other  physical  features  of  the  region  to  be  modeled.   It 

kl 
is  recommended  that  the  WES  model   be  modified  so  that  it  can  be  used 

in  conjunction  with  results  of  this  study  to  determine  actual  runup  for 

areas  where  runup  may  not  be  equal  to  elevation  at  the  shoreline  and 

where  Bretschneider ' s  method  is  not  completely  satisfactory. 
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APPENDIX  A:   NOTATION 


a  Boundary  of  region 

a  Radius  of  vertical  circular  cylinder 

A  Coefficient  in  wave  height  equation 

b(a))  Amplitude  of  frequency  component   co 

b  Incident  wave  amplitude 
o 

B  Coefficient  in  wave  height  equation 

d  Water  depth,  ft 

D  Number  of  years 

F  Frequency  per  year  of  tsunami  occurrence 

2 

g  Acceleration  due  to  gravity,  32.2  ft/sec 

h  Elevation  of  the  maximum  tsunami  wave  crest  above  mean  sea  level 
at  the  shoreline 

H  Hankel  function  of  the  first  kind  of  order  n 
n 

i  v^T 

k  Wave  number 

m  Number  of  components 

n  Integer 

n  Unit  normal  vector  outward  from  region  R 
a 

P  Probability 

r  Spherical  coordinate,  ft 

R  Region  containing  the  islands 

R  Real  number 
e 

t  Time,  sec 

x  Cartesian  coordinate,  ft 

y  Cartesian  coordinate,  ft 

a  Unknown  coefficient 
n 

3  Unknown  coefficient 
n 

n  Wave  amplitude 

0  Spherical  coordinate,  radian 

£  Response  of  harbor 

p(co)  Phase  angle 

2 

<f>  Total  velocity  potential,  ft  /sec 


Al 


2 
<j>       Total  velocity  potential  evaluated  on  boundary  a  ,  ft  /sec 
a  2 

<f>       Velocity  potential  of  incident  wave,  ft  /sec 

2 
<j)0      Far-field  velocity  potential,  ft  /sec 

a)       Angular  frequency,  radians/sec 

V       Gradient  operator ,  ft 

§       Line  integral 


A2 
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